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Abstract 

In ultrarelativistic nuclear collisions heavy quarks are produced out of thermal equilibrium in the 
very early stage of the reaction and their thermalization time was expected to be considerably 
larger than that of light quarks. On the other hand, a strongly-interacting QGP can be envisaged 
in the heavy quark sector due to the presence of heavy-light hadron-like resonances in the 
QGP for temperatures up to ~2Tc. We investigate the consequences of such states employing 
a relativistic Langevin simulation. Hadronization is modeled by a coalescence-fragmentation 
scheme. We present the predictions for the nuclear modification factor and elliptic flow of D and 
B mesons at LHC energies and compare the results with the successful predictions of this model 
for RHIC conditions. We find similar heavy-quark thermalization effects at LHC and RHIC. 
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1. Introduction 

For heavy quarks, their mass mQ^T c ~180 MeV (critical temperature). Hence their 
kinetic equilibration time is expected to be larger than that for light partons, thus show- 
ing better sensitivity to the in-medium interactions. Due to a slower thermalization a 
smaller elliptic flow of D mesons, , was expected; however, it was suggested that a siz- 
able could result from that of the light quarks [1] through a coalescence hadronization 
mechanism [1,9]. Surprisingly, data from the Relativistic Heavy-Ion Collider (RHIC) for 
single electrons (e ± ) associated with semileptonic B and D decays in semi-central Au-Au 
collisions exhibited a i>2 of up to 10% [2,5], indicating substantial collective behavior of 
charm (c) quarks consistent with the assumption of a c-quark vi similar to the one for 
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Fig. 1. (Color online) HQ pt distributions at RHIC (left) and LHC (right). Circles are the initial 
distribution, the bands show the ones after interactions in the QGP. Dashed lines are the thermal 
distributions shown for comparison. 

light quarks, apart from a pr shift due to radial- flow effects [1]. In addition, the nuclear 
suppression factor was found to be comparable to the pion one, Raa—0-3 [3,4]. Pertur- 
bative QCD (pQCD) calculations of radiative energy loss cannot explain these findings, 
even after inclusion of clastic scattering. Based on lattice QCD (1QCD) results which 
suggest resonance structures in the meson-correlation function at moderate tempera- 
tures, an effective model for heavy-light quark scattering via D and B resonances was 
suggested [6] . When implementing such a picture into a Langevin simulation the results 
for semileptonic e spectra [7] are in reasonable agreement with RHIC data [3,4,5]. 



2. Heavy-Quark Transport in a sQGP Fireball 

We employ a Fokker-Planck approach to evaluate drag and diffusion coefficients for c 
and b quarks in the QGP based on clastic scattering with light quarks via D- and -B-meson 
resonances with a width r=400-750 MeV (supplemented by perturbative interactions) [6]. 
Heavy-quark (HQ) kinetics in the QGP is treated as a relativistic Langevin process [7], 
and the medium is modeled by a spatially homogeneous elliptic thermal fireball which 
expands isentropically. 

For RHIC the details of the fireball and the parameters used can be found in [7]; 
for LHC we extrapolate to d-/V c h/d?/~1400 for central t/snn=5.5 TeV Pb-Pb collisions. 
The expansion parameters are adjusted to hydrodynamic simulations, resulting in a total 
lifetime of Tfb~6 fm/c at the end of a hadron-gas QGP mixed phase and an inclusive 
light-quark elliptic flow of (v2)=7.5%. The main change with respect to RHIC is the 
expectation of a new "phase" in which hadron-like resonances melt at higher temperature. 
In order to implement this "melting" of D- and B-mesons above Tdiss=2 T c ~360 MeV, 
a "transition factor" (1 + exp[(T — Tdi ss )/A]) _1 (A=50 MeV) has been introduced into 
the transport coefficients. 

Initial HQ spectra are computed using PYTHIA with parameters as used by the ALICE 
Collaboration [10], and c and b quarks are hadronized into D and B mesons at T c by 
coalescence with light quarks [1]; "left-over" heavy quarks are treated by 5-function 
fragmentation. 
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Fig. 2. (Color online) Predictions of relativistic Langevin simulations for heavy qurks in an sQGP for 
b=7 fm, ^/sjvjv=5.5 TeV Pb-Pb collisions: Raa (left) and V2 (right) for D and B mesons. 



3. Results and Predictions for RHIC and LHC 



In Fig. 1 the invariant pr distributions for c (red) and b (blue) quarks at RHIC (left) 
and LHC (right) are shown. It is evident that the charm-quark distributions after rescat- 
tering in the sQGP (band corresponding to the uncertainty in the resonance width) are 
close to the thermal distribution for pt<2iti c where most of the yield resides. On the 
contrary, b quarks are predicted to be off-equilibrium (even at LHC energies). 

Fig. 2 summarizes our results for HQ diffusion in a QGP in terms of the meson Raa(pt) 
and V2(pt) for b=7 fm Pb-Pb collisions at the LHC. Our most important findings are: 
(a) resonance interactions substantially increase (decrease) (Raa) compared to per- 
turbative interactions; (b) b quarks are much less affected than c quarks, reducing the 
effects in the spectra; (c) there is a strong correlation between a large vi and a small 
Raa at the quark level, which, however, is partially reversed by coalescence contributions 
which increase both V2 and Raa at the meson (and e^) level. This feature turned out 
to be important in the prediction of spectra at RHIC; (d) the predictions for LHC 
are quantitatively rather similar to our RHIC results [7] , due to a combination of harder 
initial HQ-pT spectra (see circles in Fig. 1) and a decrease in interaction strength in the 
early phases where non-perturbative resonance scattering is inoperative. Therefore we 
conclude that if at RHIC the dominant contribution to HQ interactions are hadron-like 
resonances, at LHC we should observe an Raa and v<i pattern similar to RHIC. 
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